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Abstract 

The Hamiltonian theory of zero-curvature equations with spectral parameter on 
an arbitrary compact Riemann surface is constructed. It is shown that the equations 
can be seen as commuting flows of an infinite-dimensional field generalization of the 
Hitchin system. The field analog of the elliptic Calogero-Moser system is proposed. An 
explicit parameterization of Hitchin system based on the Tyurin parameters for stable 
holomorphic vector bundles on algebraic curves is obtained. 
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1 Introduction 



The main goal of this paper is to construct a Hamiltonian theory of zero curvature equations 
on an algebraic curve introduced in |jj , and identify them as infinite-dimensional field analogs 
of the Hitchin system 0. 

The zero curvature equation 

d t L — d x M + [L, M] = 0, (1.1) 
where L(x, t, A) and M(x, t, A) are rational matrix functions of a spectral parameter A 

L = uq(x, t) + J2 u is(%i t)(\ - \i)- s i M = v (x, t) + J2 Vjk(x, t)(X - ^)- fc , (1.2) 

i,s j,k 

of degree n and m, respectively, was proposed in [|J as one the most general type of rep- 
resentation for integrable systems. Equation ( |1 . 1|) , which has to be valid identically in 
A, is equivalent to a system of (n + m + 1) matrix equations for the unknown functions 
UQ,VQ,Ui S ,Vjk- The number of the equations is less than the number of unknown functions. 
That is due to a gauge symmetry of ( |1.1|) . If g(x, t) is an arbitrary matrix function then the 
transformation 

L ^ g^- 1 + gLg-\ M i — > g t g~ x + gMg^ 1 (1.3) 

maps solutions of ( |1 . 1| ) into solutions of the same equations. The gauge transformation can 
be used to normalize L and M. For example, in the gauge uq = vq = the numbers of 
equations and unknown functions are equal. Hence, equation (|1.1|) is well-defined. 

The Riemann-Roch theorem shows that the naive direct generalization of the zero cur- 
vature equation for matrix functions that are meromorphic on an algebraic curve of genus 
g > leads to an over-determined system of equations. Indeed, the dimension of (r x r) 
matrix functions with fixed degree d divisor of poles in general position equals r 2 (d — g + 1). 
If divisors of L and M have degrees n and m, then the commutator [L, M] is of degree 
n + m. Therefore, the number of equations r 2 (n + m — g + 1) is bigger that the number 
r 2 (n + m — 2g + 1) of unknown functions modulo gauge equivalence. 

There are two ways to overcome the difficulty in defining the zero curvature equations on 
algebraic curves. The first one is based on a choice of special ansatz for L and M. On this 
way a few integrable systems were found with Lax matrices that are elliptic functions of the 
spectral parameter. The second possibility, based on a theory of high rank solutions of the 
KP equation 0, was discovered in [|IJ. It was shown that if in addition to fixed poles the 
matrix functions L and M have moving rg poles with special dependence on x and t, then 
equation ( |1 . 1| ) is a well-defined system on the space of singular parts of L and M at fixed 
poles. Recently, an algebraic construction of the zero curvature equations on an algebraic 
curve was proposed in M. 

If matrix functions L and M do not depend on x, then (|1.1| ) reduces to the Lax equation 

d t L=[M,L\. (1.4) 

A theory of the Lax equations on an algebraic curve, was briefly outlined in [0] . In the next 
section for each effective degree N > g divisor Dona smooth genus g algebraic curve T we 
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introduce a space C D of the Lax matrices, and define a hierarchy of commuting flows on it. 
The spaces of the Lax matrices associated to equivalent divisors are isomorphic. If D = K, is 
the divisor of zeros of a holomorphic differential, then the space L K is identified with an open 
set of the cotangent bundle T*(A4) of the moduli space of semistable holomorphic vector 
bundles on T, i.e. with an open set of the phase space of th Hitchin system. The commuting 
hierarchy of the Lax equations on L K are commuting flows of the Hitchin system. 

The conventional approach to a theory of the Hitchin system is based on a representation 
of T*(A4) as the Hamiltonian reduction of free infinite-dimensional system modulo infinite- 
dimensional gauge group. In the finite-gap or algebro-geometric theory of soliton equations 
involutivity of the integrals of motion does not come for granted, as in the case of the 
Hamitonian reduction. Instead, the commutativity of the hierarchy of the Lax equations 
is a starting point. It implies involutivity of the integrals, whenever the equations are 
Hamiltonian. 

The Lax matrices provide an explicit parameterization of the Hitchin system based on 
Tyurin parameters for framed stable holomorphic bundles on an algebraic curve ||. Let V 
be a stable, rank r, and degree rg holomorphic vector bundle on T. Then the dimension of 
the space of its holomorphic sections is r = dim H°(T, V). Let ci, . . . , oy be a basis of this 
space. The vectors 0^(7) are linear independent at the fiber of V over a generic point 7 G T, 
and are linearly dependent 

E^(7.) = (1.5) 
i=i 

at zeros 7 S of the corresponding section of the determinant bundle associated to V. For 
a generic V these zeros are simple, i.e. the number of distinct points 7^ is equal to rg = 
deg V, and the vectors a s = (a l s ) of linear dependence Q1.5Q are uniquely defined up to a 
multiplication. A change of the basis Oi = J2j 9ij a j corresponds to the linear transformation 
of the vectors a s , a s = g T a s . Hence, an open set M. C M. of the moduli space of vector 
bundles is parameterized by points of the factor-space 

M = M /SL r , M C S ra (r x CP"- 1 ) , (1.6) 



where SL r acts diagonally on the symmetric power of CP 7 '- 1 . In [p], [?] the parameters (7^, a s ) 
were called Tyurin parameters. Recently, the Tyurin parameterization of the Hitchin system 
for r = 2 was found [§]. 

In section 3 we show that the standard scheme to solve conventional Lax equations using 
the concept of the Baker- Akhiezer function is evenly applicable to the case of Lax equations 
on algebraic curves. We would like to emphasize that solution of the Lax equations via the 
spectral transform of the phase space to algebraic-geometric data does not use a Hamiltonian 
description of the system. Moreover, a'priori it's not clear, why the Lax equations are 
Hamiltonian. In Section 4 we clarify this problem using the approach to the Hamiltonian 
theory of soliton equations proposed in TO, |ll|]. It turns out that for D = fC the universal 



two-form which is expressed in terms of the Lax matrix and its eigenvectors coincides with 
canonical symplectic structure on the cotangent bundle T*(M). If the divisor D/c — D — K 
is effective, then the form is non-degenerate on symplectic leaves defined by a choice of the 
orbits of the adjoint action of SL r on the singular parts of L 6 CP at the punctures P m G -D^. 
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In section 5 for each degree N > g divisor DonTa commuting hierarchy of zero curvature 
equations is defined. The infinite-dimensional phase space A D of the hierarchy can be seen 
as a space of connections d x — L(x, q) along loops in A4 . We would like to emphasize that 
A D does depend on the divisor D and not simply on its equivalence class, as in the case of 
the Lax equations. If Djq is effective, then the equations of the hierarchy are Hamiltonian 
after restriction on symplectic leaves. 

The Riemann surface of the Bloch solutions of the equation 

(d x -L(x,q))ilj(x,q)=0, x G S\ q ET (1.7) 

is an analog of the spectral curves in the x-independent case. Algebro-geometric solutions 
of the hierarchy are constructed in the last section. Note, that they can be constructed in 
all the cases independently of whether the equations Hamiltonian or not. 

It is instructive to present two examples of the zero curvature equations. The first one is 
a field analog of the elliptic Calogero-Moser system. The elliptic CM system is a system of 
r particles with coordinates g« on an elliptic curve with the Hamiltonian 

where p(q) is the Weierstrass function. In fl2j the elliptic CM system was identified with a 



particular case of the Hitchin system on an elliptic curve with a puncture. In section 5 we 
show that the zero curvature equation on an elliptic curve with a puncture is equivalent to 
the Hamiltonian system which can be seen as the field analog of the elliptic CM system. For 
r = 2 this system is equiavalent to the system on a space of periodic functions p(x),q(x) 
with canonical Poisson brackets 



{p(x),q(y)} = 8(x -y). (1.9) 



The Hamiltonian is 



f 

lx. 

2(1 -i 

The second example is the Krichever-Novikov equation 



H = f (p 2 (l - ql) ~ Q ™ + 2(1 - 3ql)p(2q)) dx. (1.10) 



13 1 

q t = -^xxx + g— (l - qlx) - ^Q(q)qI, (i-ii) 



where 

^2 



Q(q)=d q <S> + <!> 2 , ® = $(q,y)=((q-y) + ((q + y)-((2q). (1.12) 

Note, that Q(q) does not depend on y. Each solution q = q(x,t) of ( |1.11|) defines a rank 2, 
genus 1 solution of the KP equation by the formula 

8u(x, y, t) = (q 2 xx - l) q~ 2 - 2q xxx q- 1 + 8 fc $ + 4q 2 x - $) . (1.13) 



Equation ( 1.12| ) has zero curvature representation on the elliptic curve with puncture with 



r = 2. The difference between the two examples is in the choice of orbits at the puncture. 
In the first example the orbit is that of the diagonal matrix diag(l, —1), while the second 
example corresponds to the orbit of the Jordan cell. 



4 



2 The Lax equations 



We define first the space of Lax matrices associated with a generic effective divisor D on T, 
and a point (7, a) = {■y s , a s } of the symmetric product X = S rg (T x CP r ^ 1 ). Throughout 
the paper it is assumed that the points 7 S G T are distinct, 7s 7^ 7fc. 

Let T^^ a be the space of meromorphic vector functions / 011 T, that are holomorphic 
except at the points 7 S , at which they have a simple pole of the form 

M = ^y+0(D, (2.1) 
The Riemann-Roch theorem implies that 

dim J- Jt a > r(rg — g + 1) — rg(r — 1) = r. (2-2) 



The first term in ( |2.2| ) is dimension of the space of meromorphic vector-functions with simple 
poles at 7 S . The second term is the number of equations equivalent to the constraint that 
poles of / are proportional to the vectors a s . 

The space JF S of meromorphic functions in the neighborhood of 7 S that have simple pole 
at 7 S of the form (|2.1| ) is the space of local sections of the vector bundle V Jja corresponding 
to (7, a) under the inverse to the Tyurin map described in terms of Hecke modification of 
the trivial bundle. The space of global holomorphic sections of Q is just the space !F ljOL . 
Let Jid' Q be an open set of the parameters (7, a) such that dim JF 7iQ , = r. 

Let D — J2i iriiPi be an effective divisor on T that does not intersect with 7. Then we 
define a space MP a of meromorphic matrix functions M = M(q), q e T, such that: 

1°. M is holomorphic except at the points 7 S , where it has at most simple poles, and at 
the points Pi of D, where it has poles of degree not greater than mf, 

2°. the coefficient M sQ of the Laurent expansion of M at 7 S 

M= -^ + M sl + M s2 {z-z s ) +0{(z-z s ) 2 ), z s = z( ls ), (2.3) 
z z s 

is a rank 1 matrix of the form 

Ms = /xsaj <— M% = fj^al (2.4) 



where fi s is a vector. The constraint fl2.4p does not depend on a choice of local coordinate z 
in the neighborhood of j a . 

If (7, a) G then the constraints (2.4) are linear independent and 



dim M^ a = r 2 (N + rg-g + l)- r 2 g(r - 1) = r 2 (N + 1) , N = degD. (2.5) 
Central to all our further constructions is a map 

V : A/? Q .— > T 7 , a (M' ) (2.6) 
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from Af® a to the tangent space to Jid' at the point (7, a). The tangent vector d m = T>(M) 
is defined by derivatives of the coordinates 

d m z s = -tr M s0 = -a^/i s , z s = z(<y s ), (2.7) 
d m a T s = -a T s M sl + K s a T s , (2.8) 

where k s is a scalar. The tangent space to CP r ~ x at a point represented by the vector a s is 
a space of r-dimensional vectors v modulo equivalence v' — v + n s a s . Therefore, the right 
hand side of ( |2.8|) is a well-defined tangent vector to CP r ~ l . 

Simple dimension counting shows that on an open set of Ai' the linear map T> is an 
injection for N < g — 1, and is an isomorphism for iV = g — 1. Let us define the space £® a 
of the Lax matrices as the kernel of T>. In other words: a matrix function L(q) G NP a is a 
Lax matrix if 

(i) the singular term of the expansion 



z- z s 

is traceless 



+ L sl + L s2 (z - z s ) + 0((z - z 8 ) ), L s0 =(3 s a s , z s = z(>y s ), (2.9) 



a T 8 P a = tr L s0 = 0; (2.10) 

(ii) is a left eigenvector of the matrix L s \ 

a s r L sl = o^k s . (2.11) 

For a non-special degree N > g divisor D and a generic set of the parameters (7,0), the 
space C® a is of dimension 

dim C® a = r 2 (N + l)-rg- rg(r - 1) = r 2 (N - g + 1) . (2.12) 

A key characterization of constraints( [2.9| - f2".ll| ) is as follows. 

Lemma 2.1 A meromorphic matrix-function L in the neighborhood U of ^s with a pole at 
7 S satisfies the constraints (\2.16 ) and \2. 1 1\ ) if and only if it has the form 

L = ^ s (z)L s (z)^; 1 (z), (2.13) 

where L s and Q s are holomorphic in U , and det $ s has at most simple zero at 7^ . 

Proof. Let g s be a constant non-degenerate matrix such that 

a T s g s = e T x , ef = (1,0,0, ... ,0). (2.14) 

If L satisfies ( P79| , |2.10p , then, the coefficient L' s0 of the Laurent expansion at 7 S of the gauge 
equivalent Lax matrix 

L' 8 = g; 1 Lg 8 = ^- + L' 8l + 0{z-z a ), z s = z( ls ), (2.15) 

Z - Zs 
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equals fej, where / = g s 1 /3 S . Therefore, it has non-zero entries at the first column, only. 

(L' s0 r=0, j = 2,...,r. (2.16) 

Further, the vector ej is a left eigenvector for L' sl corresponding to the eigenvalue k s . Hence, 
the first row of L' sl equals 

(L' al ) u = K., (I'/ = 0, J = 2,...,r. (2.17) 
From (|2.16|J2.17| ) it follows that the matrix L s = f~ 1 L' s f s , where f s is the diagonal matrix 

/ s (*) = diag{(z-z s ),l, !,...,!}, (2.18) 



is regular at j s . Hence, the Lax matrix L has the form (|2.13|) , where 

*. = 9.f.(z). (2.19) 

Conversely suppose L has the form ( |2.13| ), and let a s be the unique (up to multiplication) 
vector such that aj$> s (z s ) = 0. Then the Laurent expansion of L at 7 S has the form ( |2.9| ). 
The trace of L is holomorphic, which implies ( 2.10|) . Using the equality aj$> s (z s ) = we 
obtain that ajL is holomorphic at 7 S and its evaluation at this point is proportional to aj. 
This implies fl2.11|) and the Lemma is proved. 

Let [D] be the equivalence class of a degree N > g divisor D. Then for any set (7, a) 
there is a divisor D' equivalent to D that does not intersect with 7. Constraints ( |2.10|) and 
( P- 11| ) are invariant under the transformation L — > hL, where h is a function holomorphic in 
the neighborhood of 7^. Therefore, the spaces C® a and C® a of Lax matrices corresponding 
to equivalent divisors D and D' are isomorphic. They can be regarded as charts of a total 
space £' D 1, the Lax matrices corresponding to [D]. 

Let us consider in greater detail the case D = /C, where /C is the zero divisor of a 
holomorphic differential dz. Then Ldz, where L e £^ a , is a matrix valued one- form that is 
holomorphic everywhere except at the points j s . The constraints (|2.10| , p.llf ) imply that the 
space T s of local sections of V 1)0 is invariant under the adjoint action of L, 

f E J- s 1 > L T (z)f(z) G F s . (2.20) 

Therefore, the gauge equivalence class of the matrix valued differential Ldz can be seen as 
a global section of the bundle End(V ltOC ) <S> fl l,0 (T). It is basic in the Hitchin system theory, 
that the space of such sections, called Higgs fields, is identified with the cotangent bundle 
T*(M). 

It is instructive to establish directly the equivalence 

C K /SL r = T*{M), (2.21) 

using the map ( |2.6j ). The formula 

(L, M)=~Y, res 7s Tr (LM) dz (2.22) 
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defines a natural pairing between CJ^ and MF a - For a generic degree (g — 1) divisor Z) the 
map (|2.6|) is an isomorphism. Therefore each tangent vector w = (z s , a s ) to Ai' at the point 
(z s = zi^fs), a s ) can be represented in the form T>(M). From ( |2.7| , |2.8| ) it follows that Q2.22| ) 
actually defines a pairing between C£ a and the tangent space T 7iQ (.Mo) 



(L, w) = J2( K sZ s + atsPs)- 



(2.23) 



This formula shows that the vector f} s and the eigenvalue k s in ( [2.10| , 2.11]) can be regarded 
as coordinates of a cotangent vector to A4' . Note, that k s under the change of dz to another 
holomorphic differential dz\ get transformed to k' s = K s dzjdz\. Therefore, the pair (7^, k s ) 
can be seen as a point of the cotangent bundle T*(T) to the curve T. 

The pairing ( |2.23| ) descents to pairing of C K /SL r with tangent vectors to M.. Indeed, 
tangent vectors to M. at a point represented by gauge equivalence class of a are identified 
with vectors a s modulo transformation — > aj + a^W, where If is a matrix. Under this 
transformation the right hand side of (|2.23 ) does not changes due to the equation 



rg 



^2P s a s = res 7s Ldz = 0, 



(2.24) 



s=l 



which is valid, because Ldz is holomorphic except at 7. 

The induced pairing of C K /SL r with T(A4) is non-degenerate. Indeed, if w = T>(M), 



then ( 2.22 ) implies that 



(L, w) = (L, M) =Y, respTr (LM) dz. 



(2.25) 



Therefore, if ( |2.23| ) is degenerate then there is a nontrivial L which has zero of order m ; at 
all the points Pi of D. That is impossible because D is a generic degree (g — 1) divisor. 

Our next goal is to introduce an explicit parameterization of C K . Recall, that we always 
assume (7, a) G Ai' . 

Lemma 2.2 The map 

L e C K 1— > K,/? s , 7s ,k s }, (2.26) 
where pairs of orthogonal vectors (a^/3 s ) = are considered modulo gauge transformations 

a s ^X s a s , Ps^K'Ps, (2.27) 



and satisfy equation ( )2.24) , is one-to-one correspondence. 



Proof. Suppose that images of L and V under ( p.26| ) coincide, then (L — L')dz is a holomor- 
phic matrix valued differential (p such that 

ot T Mls) = 0. (2.28) 

Let T^ a be the space of meromorphic vector functions with poles at 7 S of the form (|2.1| ) 
and with simple pole at a point P 6 T. By the definition of A4' , the constraints ( p.l|) are 
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linearly independent. Therefore, has dimension 2r, and the vectors of singular part 



7, a 

of / G J 7 ^ at P span the whole space C r . From (|2.28|) it follows that if / G .7\f a , then 
the differential f T ip has no poles at j s . As the sum of all the residues of a meromorphic 
differential equals zero, then the f T ip is regular at P. That implies (p(P) = 0. The point P 
is arbitrary, therefore ( 2.26|) is an injection. 



The map ( |2.26| ) is linear on fibers over (7, a). Therefore, in order to complete a proof 
of the lemma, it is enough to show that dimension of C^ a is greater than or equal to 
the dimension d of the corresponding data (/3 s ,k s ). The vectors f3 s are orthogonal to a s . 
Therefore, d equals r 2 g minus the rank of the system of equations ( |2.24j) . 

Let us show, that if (7, a) G A4' then the vectors a s span C. Suppose that they span 
an /-dimensional subspace, then by a gauge transformation we can reduce the problem to 
the case when the vectors a s have the (r — I) vanishing coordinates, a\ = 0, i > I. The 
Riemann-Roch theorem then implies, that dimension of the corresponding space J-" 7jQ , is not 
less than l(rg — g + 1) — rg(l — 1) + (r — /) = (r — l)g + r. 

If the rank of a l s is r, then equations Q2.24Q are linearly independent by themselves, but 
one of them is already satisfied due to the orthogonality condition for /3 S , which implies 
Tr (/3«aJ) = 0. Therefore the dimension of the fiber of data ( |2.26| ) over (7, a) G A4' equals 

On the other hand, for L G C^ a among constraints ( |2.10| ) there are at most (rg — 1) 
linearly independent, because a meromorphic differential can not have a single simple pole. 
Hence, dimension counting as in fl2.5| ) implies dim C^ a > r 2 (g — 1) + 1 and the Lemma is 
proved. 

Example. Let T be a hyperelliptic curve defined by the equation 

2g 

y 2 = R( X ) = X 2 9+l + J2 UiX\ (2.29) 

i=0 

A set of points 7 S on T is a set of pairs (y s , x s ), such that 

y 2 = R(x s ) . (2.30) 
A meromorphic differential on F with residues (/3 s aJ) at 7 S has the form 



dx 



'g-l rg 



i=0 s =l 




where Lj is a set of arbitrary matrices. The constraints ( 2.11|) are a system of linear equations 
defining Lf 

E a l L i x l + J2( a nPs)a? Vn Vs = w^off, n = l,...,rp. (2.32) 

i=0 s^n x " Xs 



in terms of data ( |2.26| ). In a similar way the Lax matrices can be explicitly written for any 
algebraic curve using the Riemann theta-functions. 
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For g > 1, the correspondence ( |2.26[ ) descends to a system of local coordinates on C K J SL r 
over an open set A4q of -M' , which we define as follows. 

As shown above, for (7, a) G A4' the matrix a l s is of rank r. We call (7, a) a non-special 
set of the Tyurin parameters if additionally they satisfy the constraint: there is a subset of 
(r + 1) indices S\, . . . , s r+ i such that all minors of (r + 1) x r matrix oi\, are non-degenerate. 
The action of the gauge group on the space of non-special sets of the Tyurin parameters Ai 
is free. 

Let us define charts of coordinates on a smooth bundle of equivalence classes of Lax 
matrices over M.$. Consider the open set of M.q such that the vectors ctj, j = 1, . . . , r, 
are linearly independent and all the coefficients of an expansion of a r+ i in this basis do not 
vanish 

r 

Ew, cj^O. (2.33) 



a 



r+l 



s=l 



Then for each point of this open set there exists a unique matrix W G GL r , such that aJW 
is proportional to the basis vector ej with the coordinates e* = <5], and aJ +l W is proportional 
to the vector eo = J2j e j- Using the global gauge transformation defined by W 

B s = W~ 1 f3 s , A s = W T a S} (2.34) 

and the part of local transformations 

A s -> \ S A S - B s -> A; 1 ^, (2.35) 

for s = l,...,r + l, we obtain that on the open set of Aio each equivalence class has 
representation of the form (A S ,B S ) such that 

Ai = a, i = 1, . . . , r; A r+1 = e . (2.36) 

This representation is unique up to local transformations ( ^.35[ ) for s — r + 2, . . . , rg. 

In the gauge ( |2.36|) equation fl2.24|) can be easily solved for Bi, ... , B r+ i. Using (|2.36| ), 
we get 

B ) + B Ui = - E ^ • (2-37) 



The orthogonality condition of Bj to Aj = ej implies that Bj = 0. Hence, 

rg 

B Wi = - E ^4- (2-38) 



The sets of r(g — l) + l) pairs of orthogonal vectors A s , B s modulo the transformations ( |2.35| ), 
and points {7^, k s } G S rg (T*(T)) provide a parameterization of an open set of T*(M). Here 
and below M = M$/SL r . 

In the same way, taking various subsets of (r + 1) indices we obtain charts of local 
coordinates which cover T*(Ai). In section 4 we provide a similar explicit parameterization 
of C D for divisors D such that = D — K, is an effective divisor. 
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Our next goal is to construct a hierarchy of commuting flows on a total space C of a 
vector bundle over an open set of A4 . Let us identify the tangent space T L (C D ) to CP at 
the point L with the space of meromorphic matrix functions spanned by derivatives d T L\ T=Q 
of all one-parametric deformations L(q, r) e C D of L. 



Lemma 2.3 The commutator [M, L] of matrix functions L e C® a and M e NP a is a 
tangent vector to CP at L if and only if its divisor of poles outside the points 7 S is not 
greater than D. 

Proof. First of all, let us show that the tangent space T L (C D ) can be identified with a space 
of matrix functions T on V with poles of order not greater than rrii at Pi, and double poles 
at the points 7 S , where they have expansion of the form 

T = z SJ ^- + ^ + +Tsl + 0(z- z s ). (2.39) 

(z - z s y z- z s 

Here z s is a constant, and 6> s , $ s are vectors that satisfy the constraint 

a T J s + a T s (3 s = 0. (2.40) 

The vectors a s ,j3 s are defined by L. In addition it is required that the following equation 
holds: 

aPTsi = a s n s + a s k s - ajL s i - z s ajL s2 , ( 2 -41) 
where L s2 and defined by (p.9|J2.11p, and constant. 



Equations (|2.40| ) and ( p.41| ) can be easily checked for a tangent vector d T L\ T=0 , if we 
identify (z s ,a s ,$ s ) with 

z 8 = <9 T z(7 s (r))| T=0 , a s = <9 r a s (r))| T=0 , $ 8 = <9 T /3 s (r))| r=0 (2.42) 

and T s i with 

T A = (d T L sl - z s L s2 ) | T=0 (2.43) 

Direct counting of a number of the constraints shows that the space of matrix functions that 
have poles of order rrii at Pi, and satisfy ( p.39| - |2~4l"l) equals r 2 (N + 1), which is the dimension 
of C D . Therefore, these relations are necessary and sufficient conditions for T to be a tangent 
vector. 



From ( |2.10| , |2.11| ) it follows that, if we define z s and a s with the help of formulae ( 2T7 



then the expansion of [M, L] at 7 S satisfies the constraints ( p.39[ - pl4T]) . The Lemma is 
thus proved. 

The Lemma directly implies, that the Lax equation L t = [M, L] is a well-defined system 
on an open set of C D , whenever we can define M = M(L) as a function of L that outside of 
the points 7 S commutes with L up to a meromorphic function with poles at the points Pi of 
order not greater than m ; . 
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Let us fix a point P G T and local coordinates w in the neighborhoods of the punctures 
Po, Pi G D. Our next goal is to define gauge invariant functions M a (L) that satisfy the 
conditions of Lemma 2.3. They are parameterized by sets 

a = (Pi,n,m), where n > 0,m > — m« are integers. (2.44) 



As follows from ( |2.5|) , for generic L G C° there is a unique matrix function M a (q) such 



7 

that: 

(i) it has the form ( f2~3"| , |2T4]) at the points 7 S ; 

(ii) outside of the divisor 7 it has pole at the point Pi, only, where the singular part at 
M a coincides with the singular part of w~ m L n , i.e. 

M~ = M a (q) - w- m L n (q) = 0(1) is regular at P i; (2.45) 

(iii) M a is normalized by the condition M a (Po) = 0. 
Theorem 2.1 The equations 

d a L=[M a ,L], O a = d/dt a (2.46) 

define a hierarchy of commuting flows on an open set of CP , which descents to the commuting 
hierarchy on an open set of C D /SL r . 

By definition, M a only depends on L, i.e. M a = M a (L). Equation ( |2.45|) implies that [M a , L] 
satisfies the conditions of Lemma 2.3. Therefore, the right hand side of ( p.46|) is a tangent 
vector to CP at the point L. Hence, ( j2.46|) is a well-defined dynamical system on an open 
set of C D . 



The Laurent expansion of (|2.46| ) at 7 S shows that the projection K*{d a ) G T(A4 ) of the 
vector d a G T(C D ) equals 

Tc*(d a )=V{M a ). (2.47) 

Now let us prove the second statement of the theorem. Commutativity of flows ( |2.46| ) is 
equivalent to the equation 

d a M b - d b M a - [M a , M b ] = 0. (2.48) 
The left hand side of (|2.48| ) equals zero at Pq, and, as follows from (|2.47f) its expansion at 



7 S satisfies ( p.39| - p~4T]) . Therefore, it equals zero identically, if it is regular at D. This easily 
follows from standard arguments used in KP theory. If indices a and b correspond to the 
same point Pi, i.e. a = (Pi, n, m), b = (Pi, m, mi), then in the neighborhood of Pi we have 

d a M b = w-^d a L n ^ +d a M b = w- m ^[M a ,L^] + d a M b = w~ mi [M~ , L ni ] + d a M b , (2.49) 

and 

[M a ,M b ] = [w- m L n + M~,w' mi L ni +M b ] = w~ m [L n , M b ) - w- mi [L n \ M~] + 0(1) (2.50) 
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From (|2.49 , |2.50| ) it follows that the left hand side of ( |2.48| ) is regular at Pi. From the 



definition of M a , it is regular at all the other points of D as well. In a similar way we prove 
( |2.48| ) for indices a = (Pi, n,m), b = (Pj,n', m') for Pj ^ Pj. 



Let us now define an extended hierarchy of commuting flows on generic fibers of the 
evaluation map C^ a — > L(P ) = L . Note that these fibers are invariant with respect to 
( p.46| ). Additional flows are parameterized by indices 

a = (P ,m;l), m > 0, l — l,...,r. (2-51) 

Let L be a matrix with distinct eigenvalues, and let us fix a representation of L in the form 
^ Kq^/q , where K is a diagonal matrix. Then for each L e CS a , such that L(q) = L , 
there exists a unique holomorphic matrix function = ^/q, which diagonalizes L in 

the neighborhood of q, i.e. L = ^K^~ l . For each index a of the form fl2.51|) we define M a 



as the unique matrix M a e that in the neighborhood of Pq has the form 

M a = w- m ^{w)E^- l {w) + 0(w), (2.52) 
where Ei is the diagonal matrix E\ 3 = d^d^ 1 . 

Theorem 2.2 The equations 

d a L=[M a ,L], a={P ,m;l) (2.53) 



defines commuting flows on the fiber of the evaluation map CP — > L , The flows ($.53j 
commute with flows ( \2.46 ). 



The proof is almost identical to that of the previous Theorem. 



3 The Baker- Akhiezer functions 

In this section we show that standard procedure in the algebro-geometric theory of soli- 
ton equations to solve conventional Lax equations using the concept of the Baker- Akhiezer 
functions (JT3, H |) is evenly applicable to the case of Lax equations on algebraic curves. 

Let L G C D be a Lax matrix. The characteristic equation 

r 

R{k, q) = det (k - L(q)) = ¥ + £ rj (q)k r - j = (3.1) 

3=1 

defines a time-independent algebraic curve T, which is an r-fold branch cover of T. The 
following statement is a direct corollary of Lemma 2.1. 



Lemma 3.1 The coefficients rj(q) of the characteristic equation ( \3.J\ ) are holomorphic func- 
tions on T except at the points Pi of the divisor D, where they have poles of order jrrii, 
respectively. 
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For a non-special divisor D the dimension of the space S D of sets of meromorphic functions 
{rj(Q), j — 1, . . . , r} with the divisor of poles jD equals 

dimS ° = i!!tfc±l)_ r(9 _l). (3.2) 

Note, that dimension counting in the case of the special divisor K, gives 

dim S K = r 2 (g - 1) + 1. (3.3) 

Equation (|3.1|) defines a map C D i — > S* 13 . The coefficients of an expansion of rj in some 
basis of S D can be seen as functions on CP . The Lax equation implies that these functions 
are integrals of motion. Usual arguments show that they are independent. These arguments 
are based on solution of the inverse spectral problem, which reconstruct L, modulo gauge 
equivalence, from a generic set of spectral data: a smooth curve T defined by {rj} G S D , 
and a point of the Jacobian J(T), i.e. an equivalence class [7] of degree g + r — 1 divisor 7 
on T. Here g is the genus of V. 

For a generic point of S the corresponding spectral curve T is smooth. Its genus g can 
be found with the help of the Riemann-Hurwitz formula 2g — 2 = 2r(g — 1) + deg u, where v 
is the divisor on T, which is projection of the branch points of T over T. The branch points 
are zeros on T of the function dkR{k, Q). This function has poles on all the sheets of V over 
Pi of order (r — l)mj. Because the numbers of poles and zeros of a meromorphic function 
are equal then degu = Nr(r — 1) and we obtain that 

jW(r - 1) . . 

9 = ^~ 2 L + r(g-l) + l. (3.4) 

Moreover, a product of dkR on all the sheets of T is a well-defined meromorphic function 
on r. Its divisor of zeros coincides with v and the divisor of poles is r(r — 1)D. Therefore, 
these divisors are equivalent, i.e. in the Jacobian J(T) of V we have the equality 

[v] = r(r — 1)[D] G J(T). (3.5) 

For a generic point Q = (q, k) of T there is a unique eigenvector ip = ip(Q) of L 

L{q)^{Q) = ty(Q), (3.6) 

normalized by the condition that a sum of its components ipi equals 1, 

X> = 1- (3.7) 

i=l 

The coordinates of ip are rational expressions in k and the entries of L. Therefore, they 
define ip{Q) as a meromorphic vector-function on V. The degree of the divisor 7 of its poles 
can be found in the usual way. Let ^(q),q G T, be a matrix with columns ip(Q l ), where 
Q l — (<?> h(q)) are preimages of q on T 

v&(g) = {^(g 1 ),...,W)}. (3.8) 
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This matrix depends on an ordering of the roots fcj(g) of ( |3.1| ), but the function F(q) = 
det 2 ty(q) is independent of this. Therefore, F is a meromorphic function on T. Its divisor of 
poles equals 27r*(7), where 71 : F — > T is the projection. In general position, when the branch 
points of T over T are simple, the function F has simple zeros at the images of the branch 
points, and double zeros the points 7 S , because evaluations of ip at preimages of 7 S span the 
subspace orthogonal to a s . Therefore, the zero divisor of F is z^ + 27, where 7 = 71 + • • - + 7r S , 
and we obtain the equality for equivalence classes of the divisors 

2M7)] = M + 2 [7] = 2[ 7 ] + r(r - 1)D, (3.9) 

which implies 

deg7 = degu/2 + rg — g + r — 1. (3.10) 

Let \l/o be the matrix defined by (|3.8|) for q = Pq. Normalization ( |3.7|) implies that \l/o leaves 
the co-vector e = (1, . . . , 1) invariant, i.e. 

e *o = e . (3.11) 

The spectral curve F and the pole divisor 7 are invariant under the gauge transformation 
L — > \l/g 1 L\l/o ) ^ ~ ^ ^0 "V> but the matrix \l/o gets transformed to the identity \l/o = ^- Let 
F = diag(/i, . . . , f r ) be a diagonal matrix, then the gauge transformation 

L — > FLF' 1 , i[)(Q) -»• f~ 1 (Q)Fi/;, where f{Q)=j^fMQ), (3.12) 



which preserves the normalization (|3.7|) and the equality \l/o = ^> changes 7 to an equivalent 
divisor 7' of zeros of the meromorphic function f(Q). The gauge transformation of L by a 
permutation matrix corresponds to a permutation of preimages Pq G T of Pq G T, which was 
used to define \l/o- 

A matrix g with different eigenvalues has representation of the form g = ^qF, where \l/o 
satisfy (|3.11| ) and F is a diagonal matrix. That representation is unique up to conjugation 



by a permutation matrix. Therefore, the correspondence described above L — > {T, 7, ^0} 
descends to a map 

£ D /£L r ^{f,[7]}, (3.13) 

which is well-defined on an open set of C D /SL r . 

According to the Riemann-Roch theorem for each smooth genus g algebraic curve F with 
fixed points q 1 , . . . , q r , and for each nonspecial degree g + r — 1 effective divisor 7 there is a 
unique meromorphic function ipi(Q), Q G F with divisor of poles in 7, which is normalized by 
the conditions ipi{q>) = Sj. Let ip{Q) be a meromorphic vector-function with the coordinates 
ipi(Q). Note, that it satisfies (|3.7|) . 

Let T be a curve defined by equation ( p.l|) , where 77 is a generic set of meromorphic 
functions on F with divisor of poles in jD. Then for each point q G F we define a matrix 
^(q) with the help of ( ft.8|) . It depends on a choice of order of the roots ki(q) of equation 
OD but the matrix function 

L(q) = *(q)K(q)*-\q), K(q) = diag(A; 1 (g), . . . , k r (q)), (3.14) 
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is independent of the choice, and therefore, is a meromorphic matrix function on T. It has 
poles of degree at Pj G D and is holomorphic at the points of the branch divisor v. 
By reversing the arguments used for the proof of ( |3.1U| ), we get that the degree of the zero 
divisor 7 of det \1/ equals rg. In general position the zeros 7 S are simple. From Lemma 2.1 
it follows that an expansion of L at 7 S satisfies constraints ( |2.1U| , pm] ), where a s is a unique 
up to multiplication vector orthogonal to the vector-columns of ^/(7 S ). Hence, L is a Lax 
matrix-function. 

If the points Pq used for normalization of ipj are preimages of Pq G T, then L, given by 
Q3.14T ), is diagonal at q = Po, and the correspondence {T,7} — > L descends to a map 

{f,m^C D /SL r} (3.15) 

which is well-defined on an open set of the Jacobain bundle over S, where it is inverse to 

Now, let L = L(q,t) be a solution of the Lax equations (|2.46| , |2~53| ). Then the spectral 



curve T of L(q, t) is time-independent and can be regarded as a generating form of the 
integrals of the Lax equations. The divisor 7 of poles of the eigenvector if>, defined by ([TEL 



37?D does depend on t a . 

It is now standard procedure to show that [7] evolves linearly on J(T). From the Lax 
equation d a L = [M a , L] it follows that, if if) is an eigenvector of L, then (d a — M a )ip is also 
an eigenvector. Therefore, 

(d a - M a )^j(Q, t) = f a (Q, t)i>(Q, t), (3.16) 
where f a {Q,t) is a scalar meromorphic function on T. The vector-function 

f(Q,t) = <p(Q,t)ip(Q,t), <p(Q,t)=exp(- fj / ro (g,r)dr) (3.17) 

satisfies the equations 

L(q, t)$(Q, t) = k$(Q, t), (d a - M a (q, t)) $(q, t) = 0. (3.18) 

It turns out that the pole divisor j(t) of if) under the gauge transform (P.17|) gets transformed 
to a time-independent divisor 7 = 7(0) of poles of if). All the time dependence of ip(Q,t) is 
encoded in the form of its essential singularities, which it acquires at the constant poles of 
fa- 
Let L(q,t) be a solution of the hierarchy of equations (|2.46| , |2T53|) . Here and below we 
assume that only finite number of "times" t a are not equal to zero. For brevity we denote 
the variables t a corresponding to indices ( |2.44| ) and ( |2.51| ) by £(i, n ,m) an d ~k{o,m;i)i respectively. 



Commutativity of the hierarchy implies that there is a unique common gauge transform 



ifj(Q,t) = (p(Q,t)ip(Q,t) such the ip solves all the auxiliary linear equations ( |3 . 1 S 



Lemma 3.2 Let ip(Q,t), ip(Q,0) = ijj(Q,0) be the common solution of equations ( \3. 1^) 
Then 
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1°. ip is a meromorphic function on T except at the points Pi and Pq, which are preimages 
on T of the points Pi G D and Pq on T, respectively. Its divisor of poles on T outside of 
Pi, Pq is not greater that 7; 

2°. in the neighborhood of Pi the function ip has the form 

$ = &,i(w, t) exp t { „ )W - m k n ^ , (3.19) 

where £i,i(w,t) is a holomorphic vector-function, and k = ki(q) is the corresponding root of 
equation j\3. 

3°. in the neighborhood of Pq the function ip has the form 



V> = Xl(w,t)exp (X)t(0,m;0 
V n 



;i)W~ m j, (3.20) 
where xi is a holomorphic vector-function such that evaluation of its coordinates at Pq equals 

x\{p 1 q) = v 1 . 

The function ip(Q,t) is a particular case of the conventional Baker- Akhiezer functions. As 
shown in [14| , for any generic divisor 7 of degree g + r — 1 there is a unique vector function 
tp(Q,t) which satisfy all the properties 1° — 3°. It can be written explicitly in terms of the 
Riemann theta-function of the curve T. 

Theorem 3.1 Let ip(Q,t) be the Baker- Akhiezer vector function associated with a non- 
special divisor^ on V. Then there exist unique matrix functions L(q,t), M a (q,t) such that 
equations ( \3. IdJ hold. 



As a corollary we get that the Lax operator L(q,t) G CP constructed with the help of ip 
solves the whole hierarchy of the Lax equations ( |2.46| , f2~53| ). 



4 Hamiltonian approach 

As we have seen, the spectral transform which identifies the space of gauge equivalent Lax 
matrices with a total space of a Jacobian bundle over the moduli space of the spectral curves 
does not involve a Hamiltonian description of the Lax equations. Moreover, a 'priori it is not 
clear, why all the systems constructed above are Hamiltonian. In this section we show that 
the general algebraic approach to the Hamiltonian theory of the Lax equations proposed 
in H [UJ and developed in JlT| is evenly applicable to the Lax equations on the Riemann 
surfaces. 

The entries of L(q) G LP can be regarded as functions on C D . Therefore, L by itself can 
be seen as matrix-valued function and its external derivative SL as a matrix-valued one-form 
on C D . The matrix \I/ (|3.8| ) with columns formed by the canonically normalized eigenvectors 
ip(Q l ) of L can also be regarded as a matrix function on C D defined modulo permutation of 
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the columns. Hence, its differential 5^ is a matrix-valued one-form on C D . In the same way 
we consider the differential 5K of the diagonal matrix K Q3.14 ). Let us define a two-form 



fl(q) on C D with values in a space of meromorphic functions on F by the formula 

Q(q) = Tr Nf^SL A 5^ — A (JJf) . (4.1) 



This form does not depend on an order of the eigenvalues of L, and therefore, is well defined 
on C D . Fix a holomorphic differential dz on T. Then the formula 



1 

U = ~2 
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res 7a Qdz + ^ resp^fMz , (4-2) 



defines a scalar- valued two-form on £ D . 
The equation 

5L = ^5K^I~ 1 + S^K^J- 1 + ^KS^' 1 (4.3) 

implies 

Q = 2 5 (Tr (if*- 1 **)) =25 (Tr f* -1 !^*)) . (4.4) 



We would like to emphasize that though the last formula looks simplier then (|4.1|) and 
directly shows that u is a closed two-form, the original definition is more universal. As 
shown in H, [TO] , it provides symplectic structure for general soliton equations. 



Lemma 4.1 The two-form u defined by ( \4-2\ ) is invariant under gauge transformations 
defined by matrices g that preserve the co-vector e^ = (1, . . . , 1), e^g = e$. 

Proof. If g preserves eo, then the gauge transformation 

L' = g-'Lg, = g^ (4.5) 

preserves normalization ( |3.7| ) of the eigenvectors. If h = (5g)g~ l , then from (|4.4|) it follows 
that under (|4.5|) Q gets transformed to Q' = Q + F, where 

F = -2 5 (Tr (Lh)) = -2Tr (5L A h + L h A h) (4.6) 

The additional term F is a meromorphic function on V with poles at the points 7 S and Pj. 
Therefore, the sum of residues at these points of the differential Fdz equals zero and the 
Lemma is proved. 

It is necessary to emphasize that in the generic case the form u is not gauge invariant with 
respect to the whole group SL r , because it does depend on a choice of the normalization of the 
eigenvectors. A change of normalization corresponds to the transformation ty' = ^V, V = L, 
where V = V(Q) is a diagonal matrix, which might depend on Q. The corresponding 
transformation of Q has the form: 

ft' = tt + 25 (Tr (Kv)) = Q + 2Tr (5K A v) , v = 5VV~\ (4.7) 

Here we use the equation 5v = v A v = which is valid because v is diagonal. 
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Let Vq C £ be a subspace of the Lax matrices such that restriction of 5kdz to Vq is a 
holomorphic differential. This subspace is a leaf of foliation on £ D defined by the common 
level sets of the functions defined on £ D by the formulae 

T iJ: i = respi Uz - z(Pi))' J kdzj , j = 0, . . . , (m* - dk), (4.8) 

where d, is the order of zero dz at Pi (compare with the definition of the universal configu- 
ration space in 0). Note, that although the functions ( |4.8[ ) are multivalued, their common 
level sets are leaves of a well-defined foliation on £ D . 



Lemma 4.2 The two-form uj defined by j\4-3j) restricted to Vq C £ d is gauge invariant, i.e. 
it descends to a form on V D = Vq /SL r . 

Let L G £ be a Lax matrix corresponding to the zero divisor K, of a holomorphic differential 
dz, then Ldz has poles at the points 7 S , only. Therefore, Vq = £ K . 

Lemma 4.3 The two-form uj on £ defined by the formula (U-W descends to a form on 
£ K /SL r , which under the isomorphism ( \2.21 ) coincides with the canonical symplectic struc- 
ture on the cotangent bundle T*(M). 

Proof. The first statement is a direct corollary of the previous Lemma. The second one 
follows from the equality 



res^Qdz = -2 5k s A5z s + ^2 $01 A 



(4.9) 



which can be proved as follows. Let L' s be the matrix defined by the gauge transformation 
fl2.15Q ) and let Q' s be the function defined by ( |4.1|) for L = L'. Then as shown above, 

(4.10) 



(4.11) 



res^Vl'gdz = res ls £ldz + 2res 7s Tr (5g s g s l A 5L — L 5g s g s 1 A 5g s g s 1 
From (|2.9| , [2.10| ) it follows that the second term in (|4.10| ) equals 

// = -2Tr (t(3 s 5a T s + 5(5 s a T s ) A 8g s g; 1 + (3 s a T Jg s gj 1 A 5g s g; 1 
Using the equality 5a^g s + a^5g s = 0, which follows from ( |2.14| ), we get 

// = 2Tr (SP. A 5aj) = -2 (tfaj A 5f3 s ) . 



(4.12) 



The matrix L' s under the gauge transformation L = f s 1 L' s f s , where f s = f s (z) is the diagonal 
matrix ( j2.18j ), gets transformed to a holomorphic matrix. Therefore, 



= res 7 X^ + 2res 7s Tr (5f s f' 1 A 5L' S - L' s 5f s f' 1 A 5f s f ! 



(4.13) 

The last term in (|4.13| ) equals zero because f s is diagonal. From ( |2.14| - [2~48l ) it follows that 



ves^n'dz = -2res 7s Tr (df s fT 1 A 6L') = 25 z s A 5k 



(4.14) 
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Equations ( f4.10| - [lT4l ) imply ( |4.9| ). In the coordinates A s and B s ( p.34l|2.38| ) on an open set 
of T*{M) the form cu due to (|Q6| ) equals 



rg rg 

oo = J2 Sk s ASz s + J2 SB I A SA b, 9 > 1 (4.15) 

s=l s=r+l 

and the Lemma is proved. 

Let us now consider the contribution to uj from poles of Ldz at the points P m of the divisor 
D/c = D — K,. The residue of the last term in ( |4.1| ) restricted to Vq vanishes. Therefore, 

u m = ~Yes Pm Qdz = res Pm Tr (LSWf- 1 A S^' 1 ) dz . (4.16) 

If Ldz has a simple pole at P m , then its residue L m is a point of the orbit O m of the adjoint 
action of GL r , corresponding to the fixed singular part of kdz, which defines the leaf Vq . 
Let £ be a matrix, which we regard as a point of the Lie algebra £ G sl r . The formula 

d^L m =[L m ,H (4.17) 

defines a tangent vector G T Lm (O m ) to the orbit at L m . The correspondence £ — > 
is a isomorphism between sl r /sl r (L m ), and T £m (O m ). Here sl r (L m ) is a subalgebra of the 
matrices, that commute with L m . Evaluation of the form (<5\I/\I/ 1 ) at is equal to £. Hence, 
Q4.16Q restricted to V® coincides with the canonical symplectic structure on the orbit O m . 
Its evaluation on a pair of vectors £, r] is equal to 

Wm (e,r/)=Tr(L m [£,7j]). (4.18) 

If Ld^ has a multiple pole at P m , then we define L m as the equivalence class of the singular 
part of Ldz. By definition two matrix differentials L and V meromorphic in the neighborhood 
of P m are equivalent if L — L is a holomorphic differential. Let Q- be a group of the 
invertible holomorphic matrix functions in the neighborhood of P m . The transformation 
L — > gLg -1 , g G defines a representation of on the finite-dimensional space of 
singular parts of meromorphic differentials. Let O m be an orbit of this representation. 

If 7i- is the Lie algebra of then the equivalence class of the right hand side of 
( |4.17l ) for £ G TL- depends only on the equivalence class of L m . Therefore, ( 14.171) defines an 



isomorphism between the tangent space to O m at L m and Tt-/Tt-(L m ), where Tt-(L m ) is 
the subalgebra of holomorphic matrix functions £ such that [L m , £] is holomorphic at P m . 
The formula 

u m = res Pm Tr (Z m [f , (4.19) 
defines a symplectic structure on O m . 

Lemma 4.4 If = D — JC > is an effective divisor, then the map 

Li — > {z s ,K s ,a s ,p s ,L m , }, (4.20) 
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is a bijective correspondence between points of the bundle CP over A4 and sets of the data 
( \4-2Q) subject to the constraints {oi^(3 s ) = 0, and 



HI 



E&«J + Yl re s Pm L m = 0, (4.21) 



modulo gauge transformations \2.2 r i ) 



If we fix a gauge on a open set of CP by Q2.36 ), then the reconstruction formulae for 
Bi, . . . , B r+ i become 

rg _ 

Bl- + i = ~ E B l A l - E Lm, (4-22) 

s=r+2 rn 

and 

B) = -Bl +1 - £) B\A{ - J2 res Pm 1% . (4.23) 

s=r+2 rn 

If g > 1, then for > the data {z s , k s , A s , B s , L m <G O m } provide explicit coordinates on 
an open set of V D . 

Theorem 4.1 Let D be a divisor such that > ; where JC is the zero divisor of a 
holomorphic differential dz. Then the form u defined by ( \4-Q ), restricted to V® descends to 
non-degenerate closed two-form on V D : 

U = Uo+ U rn, (4-24) 



where uj$ and uo m are given by is by (\j.l9j , and ( 4-1Q) , respectively. 



The representation of the form to in terms of the Lax operator and its eigenvectors provides 
a straightforward and universal way to show that the Lax equations are Hamiltonian, and 
to construct the action-angle variables. 

By definition a vector field d t on a symplectic manifold is Hamiltonian, if the contraction 
ig t uj(X) = oo(dt,X) of the symplectic form is an exact one-form dH(X). The function H is 
the Hamiltonian corresponding to the vector field dt- 



Theorem 4.2 Let d a be the vector fields corresponding to the Lax equations (\2.4Q , 2.5fy) . 
Then the contraction of u defined by (\4-2\ ) restricted to V D equals 

i da ov = 5H a} (4.25) 

where 

H a = — iesp.Tr (w- m L n+1 )dz, a= (Pi,n,m), (4.26) 

n + 1 v ' 

H a = -res Po (w- m ki)dz, a = (P ,m;l), (4.27) 
Here k\ = ki(q) is the l-th eigenvalue of L in the neighborhood of the puncture Pq. 
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Proof. The Lax equation d a L = [M a , L], d a k = 0, and equation fl3.16| ) 

d a q = M Q # + ^F a , (4.28) 
where \I/ is the matrix of eigenvectors (|3.8|) , and F a = diag(/ a (Q 1 ), • • • , f a (Q r ), imply 



i da u = -- I ^res^Adz + ^ res Pi Mz I , (4.29) 
A \s=i P.eD J 

where A = A(q) equals 

A = Tr (tf -1 [M a , L]5V - m~ l 5L(M a m + ^F a ) - ^(M^ + *F a )5ir) . (4.30) 
Using, as before, the equality L5^ — 8^K = ^5K — 5L^, we get that 

Tr (* -1 [M a , L]«J*) = Tr (^Ma^if - M a 5l}j . (4.31) 
Using the fact that K and F are diagonal, we also obtain the equation 

Tr (* _1 <JL* F ) = Tr ($# F a ) . (4.32) 



From (gggggg ) it follows that 

i d u = respTr (6K F a ) dz + R a , (4.33) 

PieD 

where 

rg 

R a = J^ res 7s Tr (8LM a ) dz + £ res P Tr (5LM a ) dz. (4.34) 
s=l p s eD 

Note that in the first term of ( [4 . 3 3| ) a sum of residues at 7 S has been dropped because K 
and F a are holomorphic at these points. 

Consider first the case of the Lax equations ( |2.46 ). The matrix M a for a = (Pi,n, m) 
is holomorphic everywhere except at the points 7 S and Pj. Therefore, Ri^ n ,m — 0. The 
corresponding diagonal matrix Fi n>m is holomorphic at the points Pj G D, j ^ i. From 
( p.45j ) it follows that Fj |n ,m in the neighborhood of Pi has the form 

F^ m = -w~ m K n + 0(l). (4.35) 

The form SKdz restricted to V D is holomorphic in the neighborhood of Pj. Therefore, 

- resp.Tr (SKFi nm )dz = res P .Tr (V m 7TW) dz = -^res P .Tr (w' m L n+1 ) dz. (4.36) 

v / n _|_ x 1 \ / 

The matrix F a corresponding to a = (Pq, m; I) is holomorphic at the points of D. Therefore, 
the right hand side of ( |4.33| ) reduces just to R a . Because, M 0)m; ; is holomorphic except at 
the points 7 S and Po, we have in this case the equation 

Po, m; i = -res Po Tr (5LM , m; i) dz, (4.37) 

which, with the help of ( |2.52| ), implies (|4.27 ). The Theorem is therefore proved. It shows 
that the Lax equations restricted to V D are Hamiltonian whenever the restriction of u is 
non-degenerate. 
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Corollary 4.1 If is an effective divisor, then the Lax equations ( \2.46\ , ^~53j ) restricted to 
V D are Hamiltonian. The corresponding Hamiltonians ( 4-26 , 4-27 ) are in involution 



{H a ,H b } = 0. (4.38) 



The basic relation which implies all equations ( |4.38|) is involutivity of all the eigenvalues of 
the Lax matrices at different points of T, i.e. 

{*i(?)A(?i)} = 0. (4.39) 

Example. Let us consider the Lax matrices on an elliptic curve T = C / {2nuo\, 2mw2} with 
one puncture, which without loss of generality we put at z = 0. In this example we denote 
the parameters 7 S and k s by q s and p s , respectively. 

In the gauge a s = e s , e{ = 5 3 s the j-th column of the Lax matrix L l] has poles only at 
the points qj and z — 0. From ( |2.10| ) it follows that LP? is regular everywhere, i.e. it is a 
constant. Equation ( |2.11|) implies that L ji (qj) = 0, % ^ j and LP 3 = pj. An elliptic function 
with two poles and one zero fixed is uniquely defined up to a constant. It can be written in 
terms of the Weierstrass cr-function as follows 

LHz) = f '('H-aM'-»W«) l + L u = Pi . (4 . 4 o) 

a(z)o-(z - qj) a(qi - qj) a{q i ) 

Let f %0 be a rank 1 matrix f %3 = a % V . As it was mentioned above, the equations oti = Ci fix 
the gauge up to transformations by diagonal matrices. We can use these transformation to 
make a 1 = b l . The corresponding momentum is given then by the collection (a 1 ) 2 and we fix 
it to the values {a 1 ) 2 = 1. The matrix L given by ( |4.40|) with f tJ = 1 is gauge equivalent to 
the Lax matrix L with a spectral parameter for the elliptic Calogero-Moser system found in 



L ll =p l} U 1 = — qj, z), i ^ j, (4.41) 

where 

*( q ,z) = ?!j?f£e«* (4.42) 
o-{z)a{q) 

Note that L has essential singularity at z — 0, which is due to the gauge transformation by 
the diagonal matrix $ = diag($(gj, z), which removes poles of L at the points qi. 

The Hamiltonian of the elliptic CM system (|1.8|) is equal to 

Ecu = \ res Tr (V 1 ^) dz. (4.43) 

For the sequel, we would like to express H C m in terms of the first two coefficients of the 
Laurent expansion of the marked branch of the eigenvalue of L at z = 0. Indeed, expansions 
of the eigenvalues of L at z = have the form 

kx(z) = (r - l)^ 1 + fcn + k 12 z + 0(z 2 ), h{z) = -z- 1 + k n + k l2 z + 0(z 2 ), I > 1. (4.44) 

The equation 

if 1 = ^ ft = TrL = ^A; / (^). (4.45) 
i=i i=i 
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implies 



#i = £*ji, E^ = 0. (4.46) 
1=1 1=1 



From (|4.44j ) and (|4.46|) it follows that 

2H C M = 2rk 12 + j2kn. (4.47) 
i=i 

Trace of L m has the only pole at z — 0. Hence, we have the equations 



res Tr (L 2 ) = 2 ^(r - l)fc n - £ fc n J = 0, (4.48) 

res Tr(L 3 ) = 3 ( (r - 1) 2 A; 12 + (r - 1)A£ + £ - fc*)^ = 0. (4.49) 
Equations ([4.48 ) and ( 4.49|) imply 



H 1 = rk u , 2H CM = r 2 k l2 + rk 2 u . (4.50) 
Our next goal is to construct the action-angle variables for uj. 



Theorem 4.3 Let L G C D be a Lax matrix, and let ^ s be the poles of the normalized (3.7) 
eigenvector ip . Then the two-form uj defined by ( \4-Q ) is equal to 

u = 9 J2 Sk{%)A5z{%). (4.51) 

8=1 

The meaning of the right hand side of this formula is as follows. The spectral curve is 
equipped by definition with the meromorphic function k(Q). The pull back to T of the 
abelian integral z(Q) = J® dz on T is a multi-valued holomorphic function on T. The 
evaluations k(%), z(j s ) at the points % define functions on the space C D , and the wedge 
product of their external differentials is a two-form on C D . (Note, that differential Sz(j s ) of 
the multi- valued function z(g s ) is single- valued, because the periods of dz are constants). 

Proof. The proof of the formula ( |4.51| ) is very general and does not rely on any specific form 
of L. Let us present it briefly following the proof of Lemma 5.1 in JlTJ (more details can be 
found in [|i~6[|). 



Let 7^, P- be preimages on T of the points -y s G T and Pj G D. Then the form to is equal 

to 

1 



rg 



uj — — - £ I £ res iVldz + £ res p iVtdz , (4-52) 

1 j=l \s=l 3 i i J 

where Vt is a meromorphic function on V defined by the formula 

&(Q) = i>*{Q)o~L(q) A 6ip(Q) - ip*(Q)5ip(Q) A 5k, Q=(k,q)ef. (4.53) 
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The expression ip*{Q) is the dual eigenvector, which is the row-vector solution of the equation 

r(Q)L(q) = W{Q), (4-54) 

normalized by the condition 

P(QMQ) = 1- (4-55) 

Note that ip*{Q) can be identified with the only row of the matrix \l/ _1 (g) which is not 
orthogonal to the column ip{Q) of ^{q). That implies that ip*{Q) as a function on the spectral 
curve has poles at the points 7^, and at the branching points of the spectral curve. Equation 
( |4.55| ) implies that it has zeroes at the poles 7 S of ip n (Q)- These analytical properties will 
be crucial in the sequel. 

The differential Qdz is a meromorphic differential on the spectral curve T. Therefore, the 
sum of its residues at the punctures P- , 7^' is equal to the negative of the sum of the other 
residues on T. There are poles of two types. First of all, Cl has poles at the poles % of ip- 
Note that 5i/j has pole of the second order at %. Taking into account that ip* has zero & t 7s 
we obtain 

res^.ft = W5IaI>)(%) A 5z(%) + 5k(%) A 5z(%) = 25k(%) A ^(7,). (4.56) 

The last equality follows from the standard formula for variation of the eigenvalue of an 
operator, ip*5Lip = 5k. 

The second set of poles of Q is the set of branch points qi of the cover. The pole of i/j* 
at qi cancels with the zero of the differential dz, dz(qi) = 0, considered as differential on 
T. The vector-function ip is holomorphic at g^. If we take an expansion of ip in the local 
coordinate (z — z(qi)) 1 ^ 2 (in general position when the branch point is simple) and consider 
its variation we get that 

5i> = -^5z( qi )+0(l). (4.57) 

dz 

Therefore, 5ip has simple pole at qi. In the similar way we have 

dk 

5k = -—5z( qi ). (4.58) 



Equalities (^4.57|) and ( [4.58|) imply that 

5kdz 



res ?i (ijj*5L A Sip) dz = res qi 



(il)*5Ldip) A 



dk 



(4.59) 



Due to skew-symmetry of the wedge product we we may replace 5L in Q4.59| ) by (5L — 5k). 
Then, using the identities ip*(5L — 5k) = 5ip*(k — L) and (k — L)dip = (dL — dk)^, we obtain 

res qi (ip*5L A 5ip) dz = -Yes qi (5ip*ip) A 5kdz = res q .(ijj*5^) A 5kdz. (4.60) 

Note that the term with dL does not contributes to the residue, because dL(qi) = 0. The 
right hand side of (|4.60| ) cancels with a residue of the second term in the sum (|4.53| ) and the 
Theorem is proved. 
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Remark. The right hand side of ( |4.51| ) can be identified with a particular case of universal 
algebraic-geometric symplectic form proposed in || . It is defined on the generalized Jacobian 
bundles over a proper subspaces of the moduli spaces of Riemann surfaces with punctures. In 
the case of families of hyperelliptic curves that form was pioneered by Novikov and Veselov 



Let (f>k be coordinates on the Jacobian J(T) of the spectral form. The isomorphism of 
the symmetric power of the spectral curve and the Jacobian is defined by the Abel map 



Hi) = E / Mi (4-61) 

s J 

where duji is the basis of normalized holomorphic differentials on T, corresponding to a choice 
of a basis of a- and 6-cycles on T with the canonical matrix of intersections. Restricted to 
V D , the differential Skdz is holomorphic. Therefore, it can be represented as a sum of the 
basis differentials 

Skdz = E Slidwi . (4.62) 

i 

The coefficients of the sum are differentials on V D of the functions 

h = <f kdz . (4.63) 

JCLi 

From (|4.51| ) it follows that ui = 5a, where 



g+r-l ~ s g 
a= £ / 6kdz = J26I i A(l) i . (4.64) 

s=l i=l 



Corollary 4.2 The form uo restricted to V D equals 



(4.65) 



For the case, when Djc > 0, this result was obtained first in fl8l . 



It is instructive to show that (|4.65| ) directly implies that uj is non-degenerate for D/c > 0. 
First of all, ( |4.65| ) implies that the forms 5I{ are linear independent. Indeed, if they are linear 
dependent at s G S D , then there is a vector v tangent to S D at s, such that 5Ii(v) = 0. Due 
to (|4.62[) we conclude d v k = 0. It is impossible for generic s, because the equation 



d v k 



R k (k,Q) 



r-j 



= 0, 



(4.66) 



implies, then, that k(Q) satisfies algebraic equation of degree less than r, i.e. the spectral 
curve T can not be r-sheeted branch cover of T. 

Second argument needed in order to complete the proof is that dimension of the space 
Sp C S D of the spectral curves corresponding to V D equals g. The number of conditions 
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that singular parts of eigenvalues of L at the points P m G are constant along V D equals 
(rdegD/c) minus 1, due to the relation 

J2 res Pm (Tr L)dz = 0, (4.67) 

which is valid, because the singular parts of L at 7 S are traceless. From ( |3.2| ) we get 

2 dim S% = Nr(r - 1) - 2r(g - 1) + 2 = 2 g = dim V D . (4.68) 



5 The zero- curvature equations 

The main goal of this section is to present the non-stationary analog of the Lax equations 
on an algebraic curve as an infinite-dimensional Hamiltonian system. 

Let A D be a space of (r x r) matrix function L(x, q) = L(x + T, q) of the real variable x 
such that: 

1°. L(x, q) is a meromorphic function of the variable q G V with poles at D and at the 



points J s (x), where it has the form (|2.3| ), i.e. L(x, q) G A/", 



rD 



3 ( x ) a T (x) 

L(x, z) = Ps{ ) \ } + L sl (x) + 0((z - z s (x)), z s (x) = z( 7 .(x)), (5.1) 

2°. the vector T>(L(x,q)) defined by map ( |2.6|) is tangent to the loop (7(x),a(x)}, i.e. 

(9 x ^ s (x) = -a? (x) (3 s (x), d x a T s (x) = -a^(x)L sl (x) + n s (x)aj(x) , (5.2) 
where k s (x) is a scalar function. 

Remark. It is necessary to emphasize, that although the loops S 1 i — > M D / SL r are lifted to 
matrix functions L'(x, q) G Af D , x E R, such that 

L'(x + T,q) =gL'{x,q)g- 1 + d x gg-\ 1 = g{x) e GL T1 (5.3) 

without loss of generality we may consider fucntions periodic in x, because V with the 
monodromy property ( |5.3| ) is gauge equivalent to a periodic matrix function L. 

The space A® of the matrix functions, corresponding to a loop a = {^(x), a(x)} in M. Q) 
is the space of sections of finite-dimensional affine bundle over the loop, because for any 
two functions L\, L 2 G their difference is the Lax matrix, L\ — L 2 G C D . Therefore, 
for a generic divisor D the space A® is non-trivial only if degD = N > g. The functional 
dimension of A^ is equal to r 2 (N — g + 1), while the functional dimension of A D equals 
r 2 (N + l). 

Lemma 5.1 If D = JC is the zero divisor of a holomorphic differential dz, then the map 

LeA K \ — ► {a s (x), (l s {x), -y s (x), k s (x)} (5.4) 
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is a bijective correspondence of A K and the space of functions periodic in x such that 

rg 

d x z( ls (x)) = -a T s {x)(3 s {x), J2^(x)a s (xf = 0, (5.5) 

s=l 

modulo the gauge transformations 

a s (x) i — > X s (x)a s (x), (5 S \ — > X^ 1 (x)(3 s (x), k s (x) i — > k s (x) + d x In X s (x) (5.6) 

a s (x) .— > W(x) T a s (x) } (3 s (x) .— > W^aO&fc), (5.7) 

where \ s (x) is a non vanishing function periodic in x and W(x) E GL r a periodic non- 
degenerate matrix function. 



Note that from ( |5.2|) it follows that locally in the neighborhood of 7 s (x) the matrix function 
L(x, Q) G A® can be regarded as connection of the bundle V over S 1 x T along the loop 
{7(2), a(x)}. Indeed, if T is a space of local sections of this bundle, which can be identified 
with the space of meromorphic vector functions f(x, z) that have the form ( |2.1| ) in the 
neighborhood of 7 S , then 

(d x + L T (x,z))f(x,z) 6f, (5.8) 



Another characterization of the constraints (15. 21) is as follows. 



Lemma 5.2 A meromorphic matrix-function L in the neighborhood of r y s {x) with a pole at 
7 S (x) satisfies the constraints ^5.2j if and only if there exists a holomorphic matrix function 
$ s (x, z) with at most a simple zero o/det$ s at 7 S such that L is gauge equivalent 

L = ^L^; 1 + d x $ s $7 1 (5.9) 

to a holomorphic matrix function L. 



The tangent space to A D is the space of functions of x with values in the tangent space to 
the space of Lax matrices T(C D ). 

Lemma 5.3 Let L e A® and M e M^ a , x ), then the commutator [d x -L, M] = M X +[M, L] 
is a tangent vector to A D at L if and only if its divisor of poles outside of r y s {x) is not greater 
than D. 



From equations ( |5.2| ) it follows that the Laurent expansion of the matrix function T 



M x + [L, M] at the point 7 s (x) has the form ( 2.39|) , where z s and a s are given by formulae 
( p.7| , p.8| ). That proves that T is a tangent vector to L D . 

Lemma 5.3 shows that the zero-curvature equation 

L t = M x + [M,L] (5.10) 

is a well-defined system, whenever we can define M(L), such that the conditions of the 
Lemma are satisfied. Our goal is to construct the zero-curvature equations that are equivalent 
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to differential equations. That requires M(L) to be expressed in terms of L and its derivatives 
in x. 

It is instructive enough to consider the case when all the multiplicities of the points 
Pi G D equal m» = 1. Let Aq be an open set in A D such that the singular part of L e Aq 
at Pi has different eigenvalues 

L(x,q) = w^dix) u {i) {x)Cr 1 {x) + 0(1), Wi = Wi(q), w^Pi) = 0, 

u (i) = diag^Or),...,^^)) , ^(aO^i^s), fc ^ Z. (5.11) 



Lemma 5.4 Lei L(x,w) be a formal Laurent series 



CO 



L = £ (5-12) 

i=-i 

swc/i tZiat l-i(x) = C(x)u(x)C~ 1 (x), where u is a diagonal matrix, with distinct diagonal 
elements. Then there is a unique formal solution \l/o — ^o{x,w) of the equation 

(d x -L(x,w))*(x,w) = 0, (5.13) 

which has the form 

* (x,v>) = C(x) (jtis(x)w s ) e ^o h{x '' w)dx \ h = diag^, . . . , h r ), (5.14) 



\s=0 



normalized by the conditions 

^ *' J . CU) 0. (5.15) 
The coefficients £ s (x) of ( \5. l^J and the coefficients h s (x) of the Laurent series 

00 

h(x,w) = ^ h s (x)w s , h-i = u, (5.16) 

s=-l 

are differential polynomials of the matrix elements of L. 

Substitution of (|544D into (gig ) gives a system of the equations, which have the form 

hs-[u,€ a+ i] = R(to,...,&,;h Q ,...,h a -i), s = -1, 0, 1, . . . . (5.17) 

They recursivelyly determine the off-diagonal part of £ s +i, and the diagonal matrix h s as 
polynomial functions of matrix elements of k(x), % < s. 

Corollary 5.1 Let ^q be the formal solution ( 5.14 ) °f equation ( S.lty) . Then for any di- 
agonal matrix E the expression w~ m ^/oE^/Q 1 does not depend on xq, and is formally mero- 
morphic, i. e. it has the form 



w-^qEVq 1 = ™ s (x)w- s . (5.18) 

s=—m 

The coefficients m s (x) are differential polynomials on the matrix elements of the coefficients 
k(x). 
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Expression ( |5.18| ) is meromorphic and does not depend on x , because the essential singu- 
larities of the factors commute with E and so cancel each other. 

We are now in position to define matrices M a , 

a = (Pi,m;l) , m > 1, / = l,...,r, (5.19) 

which are differential polynomials on entries of L, and satisfy the conditions of Lemma 
5.3. Let \I/ (x, g) = ^/ (x,w(q)) be the formal solution of equation (|5.13|) constructed above 
for the expansion ( |5.11D of L E A® at p. Then, we define M(j >m; q(x, q) as the unique 
meromorphic matrix function, which has the form ( |2.3| , |2.4| ) at the points 7 s (x), and is 
holomorphic everywhere else except at the point Pj, where 

M {iml) {x,q) = w- m (q)^ (x,q)E l ^ 1 (x,q) + 0(1), E? = 8\V l . (5.20) 

As before, we normalize M^ m .^ by the condition M(i tm; ^(x, Po) = 0. 

It is necessary to mention, that M a , as a function of L, is defined only locally, because it 
depends on a representation of the singular part of L at Pj in the form (|5.11| ). 

Theorem 5.1 The equations 

d a L = d x M a + [M a , L], a = (P,m;/) (5.21) 
define a hierarchy of commuting flows on Aq . 



Let the coefficients of ( 5.12 ) be periodic functions of x. Then, Lemma 5.4 implies that 

V (x + T,w) = ty {x,w)e p{w \ p= h(x,w)dx. (5.22) 

Jo 

Therefore, the columns of ^/q are Bloch solutions of equation ( |5.13| ), i.e. the solutions that 
are eigenvectors of the monodromy operator. The diagonal elements of the matrix p(w) are 
the formal quasimomentum of the operator ( |5.13| ). 

Our next goal is to show that for > the zero curvature equations are Hamilto- 
nian on suitable symplectic leaves, and identify their Hamiltonians with coefficients of the 
quasimomentum matrices Pi corresponding to the expansion (|5.11| ) of L at the punctures P« 

oo 

Pi( w ) = H {i,s)W s } = diag {H^.i)}, (5.23) 

s=-l 

Let us fix a holomorphic differential dz with simple zeros, and a set of diagonal matrix 
functions v^\x). Then for a divisor D, such is effective, we define first a subspace B D 
of ^4^ by the constraints 

d x (u {i \x) -v {t \x)) = 0, (5.24) 

where are the matrices of eigenvalues ( |5.11| ) of the singular parts of L E Aq. Next we 
define a foliation of B D . The leaves Vq of the foliation are parameterized by sets of constant 
diagonal matrices c ( - m * ) with distinct diagonal elements, and are defined by the equations 

u {m \x) - v im \x) = c (m) , if dz{P m ) 0, (5.25) 
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We would like to stress the difference between the constraints ( |5.24| ) and ( |5.25| ). Equations 



( f). 24 ) imply that for all the points of the divisor D the differences (u^(x) — v^'(x)) are 



x — independent matrices. For P m G we requier additionally that the difference equals 
to the fixed matrix. 



(5.26) 



As before, we define a two-form on Vq by formula ( |4.2| ), where now 

Q(q) = Tr (j*° +T (V^SL A 6%) dx - (V 1 ^) (x ) A 6pj 

and \I/ is the matrix of the Bloch solutions of (|5.13|) , i.e. 

(d x - L(x,q))V{x,q) = 0, ^(x + T,q) = ^(x,q)e p{q) . (5.27) 

We would like to emphasize that this definition is a slight modification of the formula for 
symplectic structure for soliton equations, proposed in ||. The second term in ( |5.26| ) gives 



zero contribution in the conventional theory. It is here to remove the dependence on the 
choice of Xq in the definition as may be seen as follows. The monodromy property ( |5.27|) 
implies 

Tr (V^L A 5^) (x + T) - Tr (V^L A 5^) (x) = Tr ((V 1 ^) (x) A dp) . (5.28) 
Using the equations 5L^ = 5^ x — Ldty, we obtain 

Tr (* -1 <M) = Tr (d x . (5.29) 

Hence, the form Q does not depend on a choice of the initial point x = xq. 

The same arguments as before show that u when restricted to Vq does not depend on 
the normalization of the Bloch solutions. 

Theorem 5.2 The formula (\j.2j with Q given by ( \5.2b} ) defines a closed two-form on Vq . 
This is gauge invariant with respect to the affine gauge group GL r . 

If D >K, then the contraction of to by the vector field d a defined by ^5.21\ ) equals 

i da u = 8H a , (5.30) 

where for a = (Pi, m; I) 

#(vn;0 = -resp.Tr (w~ m E l p)dz, (5.31) 
and p is the quasi-momentum matrix. 

The proof of this theorem proceeds along identical lines to the proof of the stationary analogs 
of these results presented above. First, we show that under the gauge transformation L' = 
g^Lg — g~ l d x , ^' = g^ 1 ^ the form Q gets transformed to 



Q' = Q + Tt / (2Sh A SL - 2LSh A 5h + 5h x A 8h) dx, (5.32) 
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where Sh = bgg~ l . Note that the last term does not contribute to the residues. The first 
two terms are meromorphic on T with poles at 7 S and Pj G D, only. Therefore, a sum of 
their contributions to residues of fl'dz equals zero. Hence, u does descend to a form on 

V D = V£/GL r . (5.33) 

Using (|5.32|) for the gauge transformation (|5.9|) , where $ s depends on a point z in the 
neighborhood of j s , we obtain 

rxo+T ( r \ 

res~ /s {ldz = -2 J \5k s (x) A 5z s (x) + ^ Spl(x) A 5a l s (x) J dx. (5.34) 

From ( 5. 14] ) we obtain that if dz(Pi) = 0, then 

res P mz = Tr ( H +T (5u {i \x) A [* 5u {i) )(y)dy) dx - 5u {l) (x ) A H +T 5u {i \x)dx] . 

(5.35) 

Equations (|5.24|) imply that the restriction of 5u^ to V® is x-independent. Then, from 
( p.35| ) it follows that the points Pj G K, give zero contribution to uj. From ( |5.14| ) and ( p.25| ) 
it follows that the form <5\I / \I / ~ 1 when restricted to V® is holomorphic in the neighborhood 
of P m G Djc- Therefore, in this neighborhood 

(a*,*" 1 + |^=0(l). (5.36) 

Using this equality we obtain that on Vq the following equation holds 

res Pm ttdz = -2res Pm Tr (Ld^- 1 A SW^dx^J dz. (5.37) 

Therefore, restricted to Vq the form a; is equal to the integral over the period of ( 4.24j) . The 
proof of the equation (|5.30|) , where H a is given by (|5.31|) is almost identical to the proof of 

O)- 

Important remark. The formulae ( ^.34| ) and (|5.37|) do not directly imply that to restricted to 
V is non-degenerate, because of the constraints Q5.24D . The conventional theory of the soliton 
equations, and results of the next section provide some evidence that it is non-degenerate for 
Die > 0, although at this moment the author does not know a direct proof of that. Anyway, 
equation ( |5.30| ) show that the equations ( |5.21| ) are Hamiltonian on suitable subspaces of V D . 
Then, commutativity of flows implies 

{H a ,H b } = 0. (5.38) 

The previous results can be easily extended for the case, when the leading coefficient of the 
singular part of L at the puncture p has multiple eigenvalues. 

Lemma 5.5 Let L(x,w) be a formal Laurent series ( \5.1H\ ) such that Z_i = C(x)uC~ 1 (x), 
and u = UiS 1 ^ is a constant diagonal matrix. Then there is a unique formal solution \&o = 
^0(^5^) of the equation ( \5.1c\ ), which has the form 

tfofow) = C(x) (jtUz)w~ a ) T (x^)i T(x ,w) = 1, (5.39) 
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where 

= 5 ij ; C j (x) = 0, if m = u h s > 1. (5.40) 

and the logarithmic derivative h(x,w) ofT is a formal series with non vanishing entries only 
for indices such that Ui = Uj, i.e. 

oo 

h = d x TT- x = uuj- 1 + Y, hs(x)w s , h ij = 0, if m ^ Uj (5.41) 

s=0 

The coefficients £ 8 {x) of ( \5. 35 ) and the coefficients h s (x) of (\5.4A ) are differential polynomials 
of the matrix elements of L. 

Substitution of ( |5.39| ) in ( |5.13| ) gives a system of the equations which have the form fl5.17| ) 
They recursively determines for indices such that Ui ^ Uj and the matrix h s , as 
polynomial functions of the matrix elements of h(x), % < s. 



Corollary 5.2 Let be the formal solution ( \5. 39{ ) of equation ( \5.19j) . Then for any diag- 



onal matrix E = Ei5 %3 such that Ei = Ej, if m = Uj, the expression w~ m '^oE^Q 1 does not 
depend on xq, and is formally meromorphic. The coefficients m s (x) of its Laurent expansion 
( \5.lSj ) are differential polynomials of the entries of the coefficients k(x). 

The expression w~ m ^QE^Q 1 is meromorphic and does not depend on xq because [T, E] = 0. 

The corollary, implies that if singular parts of L at the punctures Pj have multiple eigen- 
values, then the commuting flows are parameterized by sets 

a=(P t ,m;E x ), (5.42) 

where E\ is a diagonal matrix that satisfies the condition of Corollary 5.2. The Hamiltonians 
of the corresponding equations are equal to 



H a = — resp^Tr m E\ J h(x)dx^j dz. 



(5.43) 



Example. Field analog of the elliptic CM system. 

Let us consider the zero curvature equation on the elliptic curve with one puncture. We use 
the same notation as in section 4. In the gauge a s = e s , e{ = 5 3 s , the phase space can be 
identified with the space of elliptic matrix functions such that L l i has pole at the point qj(x) 
and z = 0, only. From (|5.2j ) it follows that the residue of Jji at qj equals —qj X . Therefore, 
L- 7 - 7 = pj + qj X (((z) — ((z — qj) — ((qj). Equation (|5.2j ) implies also, that L^ l (qj) = 0, i ^ j. 

Let us assume, as in the case of the elliptic CM system, that the singular part of L at 
the puncture z = is a point of the orbit of the adjoint action corresponding to the diagonal 
matrix diag(r — 1, —1, . . . , —1). Then, taking into account the momentum map corresponding 
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to the gauge tarnsformation by diagonal matrices, we get that the non-stationary analog of 
the Lax matrix for the CM system has the form 



L n - 

The values ff are fixed to 



Pi + Qix (C(z) - ({z - ft) - C(ft)) 

ajz + qj- gj) a(z - ft)q(ft) 
a(z)a(z - ft) (j(ft - ft) <r(ft) 



fifj 



i=l 

According to ( [5.34j ), the symplectic form equals 



(5.44) 
(5.45) 

(5.46) 
(5.47) 



The commuting Hamiltonians are coefficients of the Laurent expansion at z = of the 
quasimomentum, corresponding to the only simple eigenvalue of the singular part of L at 
z — 0. To find them we look for the solution of ( |5.13| ) in the form 



i[> = C{x,z)e^ h{x '' z)dx \ 

/ oo \ oo 

C= £C«(aO** , h= £ h s (x)z s , 



(5.48) 
(5.49) 



\s=0 / s=-l 

where is the eigenvector of the singular part of L, corresponding to the eigenvalue (r — V 



i.e. 



Cf = fi, (5.50) 
and the coefficients for s > are vectors, normalized by the condition 



£/<cf = o, s >o. 



(5.51) 



i=l 



Substitution of ( 5.44 , [5 .45|) into ( |5.13|) gives a system of the equations for the coordinates C, 
of the vector C: 



dxQ + hCi = q ix Ci [({z) - ({z - ft) - C(ft)] + 

fi E f& [C(z) - C(z - q 3 ) + C(ft - ft) ~ C(ft)] 

O^i 



(5.52) 



where we use the identity 

a(z + ft - ft) a(z - ftj<r(ft) 
a(z)a(z - ft) cx(ft - ft) (j(ft) 



C( z ) ~ C(z - qj) + C(ft - Qj) ~ C(?i)- 



(5.53) 



Taking the expansion of ( |5.52|) at z — 0, we find recursively the coefficients of Q and 
densities /i s of the Hamiltonians. The first two steps follows. 
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The coeffcients at z 1 of the right and left hand sides of (|5.52 ) give 



h- x = q lx + E fj = <fo + (r - ff ) = r - 1. 



The next system of equations is 



+ fiho + (r- l)cf } = pi/i + fe C; (1 } (/;Cf + /jV< 



where 



V U = C(Qj) + C(?»i) - C(?i)> ?v = ft _ 



Using (|5.51| ), we get 



(5.54) 

(5.55) 

(5.56) 
(5.57) 



Multiplying (|5.57|) by fa and taking a sum over i, we find upon using ( |5.51|) and skew- 
symmetry of V^, 

r/io = Ep</i 2 = EPi(l + g«). (5-58) 

t=i i=i 

(2) 

In the same way we get the system of equations for C\ 
rCf > + d x CP + hoCV + h x fi = p t C\ l) + q lx f l p(q l ) + h E /; (Cf^- + /jpfo-)) • (5-59) 

Consequently the expression for the density of the second Hamiltonian is 



r 2 hi 



E 



E 



fe/f p(ft) + (/*» + ft/o + E (fifidpVij + fttfpfa 



i)y?p(ft) + cf } ( /*. + »/< + E /<#ty 



(5.60) 



For the first line we have used the equation J2i (fi^ix + /feC£ ) = 0. From ( 5.57|) it follows 
that the second term in (|5.60| ) equals 



ii = -rhl + Y.\v1fi-fl + E(./r),./j\o + E fifffiVijVki ) • (s.ei ) 



For any triple of distinct integers i ^ j ^ k ^ i the following equation holds 
VijV k i + VjkVij + V ki V jk = -p(qi) - p(qj) - p(q k )- 



(5.62) 



In order to prove (|5.62 ), it is enough to check that the left hand side, which is symmetric 
function of all the variables qi,qj,qk, as a function of the variable q iy has double pole at 
qi = 0, and is regular at = qj. In the same way one can obtain the well-known relation 



VijVji = -p(qi) - p{qj) - p{qij). 



(5.63) 
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Equations ( |5.62| , p.63| ) imply 



E E fif?fkV«V u = - E r/?(r - /?)p( ft ) + E y? ] • (5-64) 



From (|5.56| ) it follows 

EOf iE ((./r)«,/j - ./r(./j).,) (C(%) - 2C( ft )) = 

= - g E ( 2r Qixx({qt) ~ E lijxxC(Qij) I + g E ~ ) C(fti)- ( 5 - 65 ) 
The first sum is equal to 

2 E ( Zrq ixx C{q t ) - E ?*i*xC(fci) J = ^ E ( 2r fe/ 2 p(ft) ~ E <lijxP(<lij) I + ^F, ( 5 - 66 ) 
where 

F = \ E ( 2r/?C(«) - E ^C(%)] • (5-67) 



The first terms in (|5.60| . |5l>4| . |5.65| ) cancel each other. The function F, as a function of the 



variable qi, has poles at the points 0, qj, j ^ i and the sum of its residues at these points 
equals 

'•/r-E'/o (5-68) 

Therefore, it has the same monodromy properties with respect to all the variables. The 
functions qi(x) represent loops on the elliptic curve. Therefore, qi(x + T) = qi(x) +b iy where 
hi is a period of the elliptic curve. The constraint ( |5.46| ) implies J2i h = 0. Then, from Q5.68| ) 



it follows that F is a periodic function of x. The densities of the Hamiltonians are defined up 
to a total derivative of periodic functions in x. Hence, a density of the second Hamiltonian 
of the hierachy equals 

r%= - + +E(^(l + fe)-4p^j) (5-69) 

- ^J2((^xxqjx- qjxxQix) ((qij) (5.70) 

+ 2 S ft 1 + 9ix){l + V*) 2 + (1 + + fe) 2 - <&) PM- (5-71) 

The transformation Pi Pi + f(x) does not change h s for s > 0. In particular, the first two 
terms in Q5.69| ) can be rewritten as 

- - (EfeC 1 + &>)) + E^ 1 + &) = ^ Efe -Pj) 2 (! + + fc*)- (5-72) 
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The symplectic form ( |5.47| ) restricted to the subspace 



(5.73) 



is non-degenerate. The Hamiltonians H s = Jq h s (x)dx restricted to this space generate a 
hierarchy of commuting flows, which we regard as field analog of the elliptic CM system. For 
r = 2 the Hamiltonian 2H\ has the form Q1.10 ), where q = q± = —q2, p — Pi — ~P2- 



6 The algebro-geometric solutions 



So far, our consideration of the Bloch solutions ( 5.27|) has been purely local and formal. For 



generic L G Aq the series ( |5.14| , |5.16j ) for the formal solutions ^(x, q), and quasimomentum 



have zero radius of convergence. The main goal of this section is to construct algebro- 
geometric solutions of the zero curvature equations, for which these series do converge and, 
moreover, have meromorphic continuations on a compact Riemann surface. 

Let T(q) be a restriction of the monodromy operator f(x) — > f(x + T) to the space of 
solutions of the equation (d x — L(x,q)f = 0, where / is a vector function. Then, we define 
the Riemann surface T of the Bloch solutions by the characteristic equation 

r 

R{fi, q) = det (// - f(qj) = // + ^ Rj{q)ff~ j = 0. (6.1) 



Lemma 6.1 The coefficients Rj{q) of the characteristic equation j\6. 1\) are holomorphic 
functions on T except at the points Pi of the divisor D. 

Proof. In the basis defined by columns of the fundamental matrix of solutions to the equation 
(d x — L)F(x, q; x ) = 0, F(x , q, x ) = 1, the operator T(q) can be identified with the matrix 

f(q) = F(x + T,q;x ). (6.2) 

A 'priory this matrix is holomorphic on T except at the points of the divisor D and at 
points of the loops 7 s (x), where L has singularities. From Lemma 5.2 it follows that in the 
neighborhood of the loop we have 

F = $ s (x, q)F(x, q)$; l (x, q), F(x , q) = 1, (6.3) 



where F is a holomorphic matrix function, and $ s is defined by ( p.l4j , ^7T8| , |2.19|) . The function 



$ s is periodic, because 7 S , a s are periodic. In the neighborhood of the loop 7 S , the functions 
Rj{q) coincide with the coefficients of the characteristic equation for F. Therefore, they are 
holomorphic in that neighborhood. The Lemma is then proved. 

It is standard in the conventional spectral theory of periodic linear operators that for 
a generic operator the Riemann surface of the Bloch functions is smooth and has infinite 
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genus. For algebro-geometric or finite-gap operators the corresponding Riemann surface is 
singular, and is birational equivalent to a smooth algebraic curve. 

It is instructive to consider first, as an example of such operators, the case, when L does 
not depend on x, i.e. L G C D . In this case the equation (d x — L)ip = can be easily solved. 
The Bloch solutions have the form 

^ = ijj e kx , (6.4) 

where ipo is an eigenvector of L, and k is the correspoding eigenvalue. These solutions are 
parameterized by points Q of the spectral curve To of L. The image of Tq under the map 
into C 1 x T defined by formula 

{k, q) G f o 1 — ► (M = e kT , q) G C 1 x T (6.5) 



is the Riemann surface T defined by ( |6.1| ), where the coefficients are symmetric polynomials 

of e *i(9)r 

For example, if r is defined by the equation 

k 2 + u(q) = 0, (6.6) 

where u(q) is a meromorphic function with double poles at the points of D, then T is defined 
by the equation 

[i 2 + 2i? ly u + 1 = 0, R^q) =cosh (Ju^j). (6.7) 



The Riemann surface defined by Q3.7D is singular. Projections onto T of the points of self- 
intersection of T are roots of the equation 

u{q) = {w) 2 > (6 - 8) 

where N is an integer. The coefficient u(q) has poles of the second order at D, u = afw~ 2 + 
0(w~ 1 ), where w is a local coordinate at Pi G D. Therefore, as \N\ — > oo, the roots of (|6.7p 
tend to the points of D. The coordinates of the singular points q^, that tend to Pj equal 

w(q ijN ) = 2Ta t (irNy 1 + 0(N- 2 ). (6.9) 

As usual in perturbation theory, for generic L each double eigenvalue qi^ n splits into two 
smooth branch points qf n . By analogy with the conventional theory we expect, that if L is 
an analytic function of x, then the differences {w^n) —w(qf N )\ < 0(N~ h ) will decay faster 
that any power of AT -1 . 



Localization of the branch points is a key element in the construction |19[ of a theory 
of theta-functions for infinite genus hyperelliptic curves of the Bloch solutions for periodic 
Sturm-Liouville operators. In p0[ a general approach for the construction of Riemann sur- 



faces of the Bloch functions was proposed. The model of the spectral curves developed in J2U 



was chosen in [21] as a starting point of the theory of general (non-hyperelliptic) infinite- 
genus Riemann surfaces. It was shown that for such surfaces many classical theorems of 
algebraic geometry take place. 
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Algebro-geometric or finite-gap operators can be seen, as operators for which there are 
only a finite number of multiple eigenvalues that split into smooth branch points. Let V 
be a smooth genus g algebraic curve that is an r-branch cover of T. Note that, unlike the 
stationary case, for given a rank r their is no relation between g, and the genus g of T. As g 
increases the dimension of the space of r-sheeted cover increases. It equals 2{g — rg + r — 1). 

Assume that the preimages P\, Pq on T of the points of a divisor D, and a point P on T 
are not branch points. The definition of the Baker- Akhiezer function corresponding to this 
data and to a non-special degree g + r — 1 divisor 7 on T is as follows: 

1°. ip is a meromorphic vector function on F except at the points P\. Its divisor of poles 
on T outside of P\ is not greater that 7; 

2°. in the neighborhood of P\ the vector function if) has the form 



= 6,1(9,*) ex P (£*(v*;Q 
\ m 



w~ m , (6.10) 



where £i,i(q,t) is a holomorphic vector-function; 

3°. evaluation of if) at the punctures Pq are vectors with coordinates (ip(Po))® = 5 d . 

Theorem 6.1 Let if)(q, t) be the Baker- Alhiezer vector function associated with a non-special 
divisor^ on Y. Then, there exist unique matrix functions M^ im .^(q,t) G MB t \ a r t \ such that 
the equations 

(% m;0 -M (iim;0 )v(g,*) = o (6.11) 

hold. 

Now, let Vi\x) be a set of periodic functions, J T v\ % 'dx = 0, and uf 1 be a set of constants. 
Then the change of the independent variables 

t(i,i;i) = xu?' ) +v®(x)+t{ i;1 . l) (6.12) 

define the Baker- Akhiezer function if), as a function of (q, t) and the variable x. From ( |6.11| ) 
it follows that 

(d x -L)if) = 0, L = Y< { u ? + fa®) M d,i;i) ( 6 - 13 ) 

i,l 

As follows from Lemma 5.2, the vector T>(M a ), a = (i,m;l), corresponding to M a under 
(|2.6p, is tangent (7 s (t a ), a s (t a )). Therefore, T>(L) is tangent to ( r y s (x),a s (a)). 

In general, L constructed above is not a periodic function of x. It is periodic, if we impose 
additional constraints on the set of data that are the curve T and the constants uf' '. We 
call the set {T,uf} admissible if there exists a meromorphic differential dp on T which has 
second order poles at P\ 

dp = -uf ] dw (w~ 2 + 0(1)) , (6.14) 
and such that all periods of dp are multiples of I'nijT 

dp=^^, m c eZ, c G #i(r, Z). (6.15) 
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Lemma 6.2 The Baker-Akhiezer function if), associated with an admissible set of data 
{r,M^' ) } satisfies the equation 

iJ;(x + T,q)=g^(x,q)fi(q), fi = e^ T , (6.16) 

where g is the diagonal matrix g = diag (p(Pq), . . . , /z(P r )). 



From ( |6.15 ) it follows that the function \x defined by the multi-valued abelian integral p is 



single-valued. Equation ( |6.16| ) follows from the uniqueness of the Baker-Akhiezer function, 
because the left and the right hand sides have the same analytic properties. 

The matrix function L constructed with the help of ip satisfies the monodromy property 

L(x + T,q)=gL(x,q)g- 1 . (6.17) 

Let S = S(T,p) be a space of curves T with meromorphic differential dp satisfying ( |6.15| ). 
We would like to mention, that the closure of S, as T — > oo, coincides with the space of all 
genus g branching covers of V. 

Corollary 6.1 A set of data T £ S, [7] £ J{f), and a set of periodic functions v^\x) define 
with the help of the corresponding Baker-Akhiezer function a solution of the hierarchy ( \5. 2\ ) 
on B D /GL r . 

The finite-gap or algebro-geometric solutions are singled out by the constraint that there is 
a Lax matrix L\ £ C nD such that 

[d x -L 1 L l }=0. (6.18) 

Indeed, let k be a function on T with divisor of poles nD, where D is the preimage of D. If 
n is big enough this exists. Let ip be the Baker-Akhiezer function on T, then as it was shown 
above there is a unique Lax matrix L\ such that 

L 1 (t,q)1>(t,q) = k(q)il>(t,q). (6.19) 

Equation ( |6.19|) implies that the spectral curve of L\ is birationally equivalent to the Riemann 
surface T of Bloch solutions for L. 

Theorem 6.2 The form uj defined by ( U-Q) and ( \5. 26]) restricted to the space of algebro- 
geometric solutions, corresponding to a set of function v\ l \x) equals 

g+r-l 

uj= £ 8p(%) A5z(%). (6.20) 



The meaning of the right hand side of this formula is analogous to that of formula ( 4.51|) 



It shows that the form to restricted to the space of algebro-geometric solutions is non- 
degenerate. 

It is well-known, that the finite-gap solutions of the KdV hierarchy are dense in the space 



of all periodic solutions (||22j|). As shown in [20 the finite-gap solutions are dense for the 



KP-2 equation as well. It seems quite natural to expect that the similar result is valid for the 
zero-curvature equations on an arbitrary algebraic curve, as well. In the conjectued scenario 
the infinite dimensional space B D can be identified with a direct limit of finite-dimensional 
spaces C nD , as n —>■ 00. We are going to address that problem in the near future. 
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